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Weighing Neutrinos with Galaxy Cluster Surveys 
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Large future galaxy cluster surveys, combined with cosmic microwave background observations, 
can achieve a high sensitivity to the masses of cosmologically important neutrinos. We show that 
a weak lensing selected sample of > 100, 000 clusters could tighten the current upper bound on 
the sum of masses of neutrino species by an order of magnitude, to a level of 0.03 eV. Since this 
statistical sensitivity is below the best existing lower limit on the mass of at least one neutrino 
species, a future detection is likely, provided that systematic errors can be controlled to a similar 
level. 



Recent experiments have placed both stringent lower 
and upper bounds on the masses of neutrinos. The lower 
bounds derive from neutrino oscillations experiments, 
which measure the difference in the squared masses of 
the coupled species. Atmospheric neutrino oscillations 
showed that at least one neutrino species has a mass 
> 0.05 eV lji implying that neutrinos make a non- 
negligible contribution to dark matter. Solar neutrino os- 
cillations showed smaller mass splittings (~ 0.008 eV 0) 
(see Fig. Q] for the allowed masses of individual species). 

The upper bounds come from cosmological large scale 
structure measurements, which are sensitive to the sum 
of the masses of all neutrinos. By combining CMB 
anisotropics from the Wilkinson Microwave Anisotropy 
Probe (WMAP) [fj and the clustering properties of 
galaxies in the Sloan Digital Sky Survey (SDSS), with 
no prior assumptions about the amplitude of galaxy bias, 
< 1.7 eV at 95% CL was obtained 0]. The limit 
can be tightened to 0.42 eV Q by incorporating weak 
lensing, Lya forest and Supernovae data, but at the ex- 
pense of introducing possible new systematics [34| ■ 

The above bounds have narrowed the allowed range 
for the sum of the masses ^ m v to within an order of 
magnitude. In this Letter, we discuss the prospect of 
closing the gap, utilizing large, future surveys of galaxy 
clusters. Such surveys, covering large fractions of the 
sky to impressive depths, are being planned in several 
wavelength bands. There are several advantages of such 
a cluster sample. First, clusters correspond to massive 
dark matter potential wells, whose abundance and spatial 
distribution are dictated by gravity alone, and have been 
simulated to high accuracy [7|. 

Second, galaxy clusters are highly clustered relative to 
the mass distribution. This bias can, at least principle, be 
determined from simulations, and it increases the signal- 
to-noise in power measurements by a factor of ~ 10— 100 
(for the relevant cluster masses below). Third, a galaxy 



cluster survey delivers several observables. As we will 
find below, the power spectrum P c (k) and abundance 
evolution dN/dz provide complementary constraints: the 
free-streaming of neutrinos suppresses fluctuation power 
on small scales, affecting the shape of P c (k), and also the 
growth rate of perturbations in the linear regime g(z), 
affecting dN/dz. The main limitation of using galaxy 
clusters is that the relation between cluster mass and the 
actual observables needs to be known to high precision, 
at least statistically. As we discuss below, this should 
be achievable using empirical calibrations and numerical 
simulations. 

Neutrino Signature. — Massive neutrinos cluster on 
very large scales, but free-stream out of small-scale dark 
matter potential wells and thus suppress density fluctu- 
ations on small scales. In addition, the evolution of per- 
turbations is no longer independent of scale. The late- 
time evolution of perturbations in a cosmology with cold 
dark matter (CDM), baryons and massive neutrinos, ac- 
counting for all post-recombination effects, can be ac- 
curately treated using a product of a scale-dependent 
growth function (Eq. 25 in Q) and a time-independent 
transfer function, reflecting conditions at the drag epoch. 
The latter is obtained numerically from CMBFAST 9] 
and includes acoustic baryonic features, which contain 
information comparable to the broader features at the 
matter-radiation equality and the sound horizon [lo|. 

Cluster Abundance and Distribution. - Since our 
analysis closely follows the treatment of 01 > here we only 
outline the methods and refer the reader to for de- 
tails. We use the results of numerical simulations from 
Jenkins et al. 7] for the differential comoving number 
density dN/dM(M, z) of clusters of total mass M at red- 
shift z (the fitting formula in their Eq. B3). Correlations 
in the spatial distribution of clusters are described by 
the power spectrum, P c (k), ta king into account redshift 
space distortions [lj] and bias [lj|. 
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FIG. 1: The total mass and three mass eigenvalues of neu- 
trinos as a function of the unknown smallest mass mi !. A 
normal hierarchy is assumed. In our analysis, we consider the 
scenario of three families of neutrinos with one massive and 
two massless. The upper dotted line shows the current most 
robust la limit The lower dotted lines show our predicted 
la constraints, for two different sets of weak lensing selection 
efficiency and completeness, from combining all three cluster 
surveys and Planck. 



Future Cluster Surveys. - Galaxy clusters can be 
identified and cataloged using several different obser- 
vational signatures. We consider three methods that 
use the X-ray emission of the hot cluster gas, the 
up-scattering of CMB photons in energy by this gas 
(Sunyaev-Zel'dovich effect, hereafter SZE), or the dis- 
tortion of the images of background galaxies by weak 
gravitational lensing (WL). To determine the detection 
mass limit M m i n (z) in each case, we adopt a flux lim- 
ited X-ray survey (guided by plans of the Dark Uni- 
verse Observatory, DUO Q, yielding a total of « 11, 500 
clusters in 6,150 deg 2 ), a temperature-decrement limited 
SZE survey (as planned with the South Pole Telescope, 
SPT H3, yielding w 20, 000 clusters in 4,000 deg 2 ), and a 
shear-limited WL survey (planned with the Large Syn- 
optic Survey Telescope, LSST 0, yielding » 200,000 
clusters in 18,000 deg 2 ). We have additionally imposed a 
floor for the minimum mass of 1O 14 /i _1 M in the X-ray 
and SZE surveys, since less massive halos correspond to 
small clusters or groups in which these two observables 
are subject to strong non-gravitational effects. If we im- 
posed the same minimum mass floor for LSST, it would 
yield « 143, 000 clusters. The mean redshifts are 0.36 for 
DUO, 0.58 for SPT and 0.43 for LSST (see H3, Fig. 1 
for the redshift distribution of clusters). To parameterize 
these effects for the clusters above the detection thresh- 
old, we introduce three additional non-cosmological pa- 
rameters in predicting both the X-ray flux and the SZ 
decrement, of the form f x (z) oc A X M^"(1 + z) 7x and 



f sz oc A SZ M^ SZ (1 + z) 7 ", respectively. Deviations from 
these power-law scalings are not captured by our param- 
eterization and would have to be addressed by allowing 
more general forms of cluster structure and evolution (see 
discussion below). For the WL survey, we assume that 
clusters (even at low-mass) produce shear signals with 
a well-defined mass-shear relation and do not introduce 
analogous free parameters. The main systematic error 
for WL surveys comes from false detections and incom- 
pleteness (discussed below). 

Fisher Matrix Forecasts. - We estimate the con- 
straints on cosmological parameters using the Fisher ma- 
trix formalism that has become a standard tool ^t|> a l~ 
lowing easy combination of different independent data 
sets and/or methods, by summing individual matrices. 
The Fisher matrix for the redshift distribution and power 
spectrum are described in detail in [llj . The Fisher 
matrix for the temperature and E-mode polarization 
anisotropy of the cosmic microwave background (CMB) 
is constructed as given by 0] . In this paper, the covari- 
ance matrix of WMAP is from the first year data, and 
the expected parameters of the Planck satellite [ljj are 
listed in Tabic I of [lj. 

The Fisher formalism requires a fiducial model around 
which variations are considered. For the cosmological pa- 
rameters, we adopt a spatially-flat, CDM model, domi- 
nated by a cosmological constant. The set of parameters 
included in our analysis is {ilt,h 2 , il m h 2 , £l„h 2 , Ode, 
w, n s , erg, t}, where all the symbols have their stan- 
dard meaning. The values are adopted from recent mea- 
surements by WMAP, as summarized in Table 1 of 0: 
{0.024, 0.14, 0, 0.73, -1, 1, 0.9, 0.17}. We find that choos- 
ing il„h 2 = 0.00058 - the current lower limit - would not 
change our results below. 

Our analysis takes into account the dominant sys- 
tematic errors for the three surveys. For X-ray and 
SZE surveys, these come from uncertainties in the mass- 
observable relation due to structure and evolution of clus- 
ters. As proposed in (ijj, EH> however, one can use clus- 
ter surveys not only as a source of cosmological infor- 
mation, but also to constrain the mass-observable re- 
lation, thereby making the survey self-calibrating. To 
account for this, we include the parameters of the mass- 
observable relation in the Fisher matrix analysis: A x 
(= 10" 4 - 159 ), p x (= 1.807) and j x (= 0) for X-ray survey; 
A sz (= 10 8 - 9 ), /3 SZ (= 1.68) and lsz (= 0) for SZE survey, 
where numbers in parentheses indicate the fiducial values 
for these parameters [2lj |. 

The main systematic limitation for WL surveys comes 
from false detections and incompleteness. In our analysis 
of the LSST-like WL survey, we used a constant shear 
S/N threshold to select clusters. False detections or miss- 
ing clusters result from statistical fluctuations in these 
ellipticities and from projections of physical structures 
along the line of sight. Several papers [i2,|2j| have done 
a comprehensive study of mass-selected clusters using 
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TABLE I: Estimated Constraints on Total Mass of Neutrinos 
(in unit of eV) . The DUO-like and SPT-like surveys are self- 
calibrated by including three non-cosmological parameters in 
our analysis. The errors in the LSST-like survey incorporate 
calibration with numerical simulations, assuming an efficiency 
of e = 40% and completeness c = 60%; numbers in parenthe- 
ses assume e = 75% and c = 70%. 'W and 'P' stand for 
WMAP and Planck respectively. 





DUO 


SPT 


LSST 




(6,150 deg 2 ) (4,000 dcg 2 ) 


(18,000 dcg 2 ) 


Pc{k) 


1.4 


1.1 


0.71 (0.42) 


P c (k)+dN/dz 


0.70 


0.72 


0.53 (0.32) 


P c (k)+C e (W) 


0.22 


0.20 


0.15 (0.11) 


P c (k)+C e (W)+dN/dz 


0.16 


0.15 


0.12 (0.10) 


Pc{k)+d(P) 


0.11 


0.10 


0.086 (0.061) 


P c (k)+C e {P)+dN/dz 


0.071 


0.062 


0.040 (0.027) 


DUO+SPT+LSST+Planck 


0.034 (0.025) 



N-body simulations. They point out that because the 
simulations depend only on gravity, the expected clus- 
ter distribution, including false detections and missing 
clusters, can be reliably calculated for any cosmological 
model. Thus, false detections and missing clusters can 
be accounted for, and their presence serves only to in- 
crease the statistical error. We here adopt the efficiency 
e = 40% and the completeness c = 60% for a 4.5 standard 
deviation detection threshold ji^- To account for this, 
we multiply our parameter error bars for the WL survey 
by a correction factor of v/[(l/e — 1) + l/e]/c ~ 2.6. 

Results and Discussion. - Table [I] shows the neu- 
trino mass constraints from our Fisher matrix analy- 
sis for CMB, and three types of cluster surveys, in- 
cluding self-calibration for X-ray and SZE surveys, as 
well as false detection and completeness for the WL sur- 
vey. The constraint from the forecast for Planck alone is 
o-(£m„) = 0.23 eV 0. 

For all three cluster surveys, the power spectrum, 
P c (fc), is a much more sensitive probe of the neutrino 
mass than the counts, dN/dz, as expected, and, by it- 
self, yields a constraint of cr(^2 mj ~ 1 eV. Combin- 
ing with WMAP and Planck improves this by factors of 
« 5 and w 10, respectively. Similarly, combining P c (k) 
with dN/dz yields a factor of 2 improvement for DUO 
and SPT, but only a modest change for LSST. (This is 
because LSST measures a smaller fraction of clusters at 
higher redshift, as seen from Fig. 1 of therefore mak- 
ing dN/dz less relevant for this survey.) Hence we find 
that each cluster survey, combined with WMAP, yields 
a(J2 m u) ~ 0-1 eV. When combined with Planck, each 
survey gives a(J2 m v) = 0-04 — 0.07 eV, very near the 
interesting limit of 0.05 eV from the Super-Kamiokande 
atmospheric neutrino experiment 

For the X-ray and SZE surveys, one can ask how much 
could be gained by completely eliminating the system- 
atics due to cluster structure and evolution. We can 



estimate this by dropping the corresponding parame- 
ters, A x , f3 x , etc. from the Fisher matrix. When com- 
bined with Planck, each survey would give <j{^2m v ) = 
0.02 — 0.03 eV, which is well below the Super-K bound. 
While a complete elimination of systematic errors is un- 
realistic, it is pointed out in [5^ that including the ad- 
ditional information from the shape of the mass func- 
tion (dN/dM) allows one to largely eliminate ambigui- 
ties caused by a wide range of possible evolutions of the 
mass-observable relation. Here we consider reducing the 
degradation of constraints due to self-calibration by com- 
bining DUO and SPT together (i.e. statistically; no over- 
lap of the survey areas is then needed). Once again with 
Planck, we now find <rQ2 m u) k 0.054 eV, back below 
the lower limit of m v ~ 0.058 eV implied by the com- 
bination of Super-K and solar oscillation experiments. 

For the WL survey, it has been demonstrated that 
a tomographic analysis can significantly reduce projec- 
tions. Also, has shown that red galaxies can be used 
to selected clusters optically with a < 5% false detec- 
tion rate. LSST will contain both lensing and optical 
observations of clusters, and projection effects and false 
detections can be removed rather than statistically sub- 
tracted. A detailed study of these and other possible en- 
hancements, such as applying a series (rather than a sin- 
gle) smoothing filter to the WL shear maps, which were 
not considered by j^, is needed. If a completeness of 
70% and efficiency of 75% can be attained without intro- 
ducing a new systematic bias, then a 0.027eV sensitivity 
of neutrinos from LSST will be achievable. In Table [i] 
numbers in parentheses in the LSST column correspond 
to c = 70% and e = 75%. 

We have checked the effects of several other system- 
atics, and found them to be small: possible redshift 
evolution in the bias b(z) [lfl Hll | ; scatter in the mass- 
observable relation j2^| ; baryon cooling altering the den- 
sity profile of clusters [28j. One important systematic 
error is the uncertainty of photometric redshifts. The 
redshift uncertainty of one single red galaxy is 0.03 for 
low redshifts and approximately twice worse for high red- 
shifts jl^l . On the other hand, a cluster of mass as large 
as 10 14 A -1 Mg should have at least 25 red galaxies [3fJ, 
and the redshift uncertainty of each cluster as a whole 
would be reduced after averaging the redshifts of the con- 
stituent galaxies. We find that an error of Az = 0.03 
would degrade the constraint on from P c (k) alone 

by a factor of ~ 2, as modes with a relatively large fcy 
will be swamped by shot noise and therefore give no 
leverage jU- However, after combining with dN/dz and 
CMB, the effect is much smaller: a « 20% increase on 
the final errors. 

There are other cosmological parameters that could 
be added in our analysis, for example, running of the 
scalar spectral index (dn s /dk) or an evolution of the 
dark energy equation of state (dw/dz). However, the 
constraint on neutrino mass derives from the evolving 
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scale-dependence of the growth function; an unevolving 
curvature in P c (k) as a function of scale, or a smooth 
evolution of the equation of state, can not mimic such 
a scale-dependent evolution. As an example, we have 
explicitly verified that adding an additional free param- 
eter, a = dn s /d\ak increases our final error on the 
neutrino mass by less than 2 percent. 

To conclude, taking into account self-calibration for 
X-ray and SZE surveys, as well as completeness and effi- 
ciency for WL surveys, we find a combined constraint for 
DUO + SPT + LSST + Planck of a(£ rn v ) = 0.034 eV. 
This implies at least a 1.7a detection of neutrino dark 
matter; improving WL selection efficiency and complete- 
ness would increase the significance to 2.3c (see Table 
I)- 

We consider how our constraints could be improved 
by a measurement of the power spectrum from the 
Lya forest. Using the one-dimensional power spectrum 
k\\Pm{k\\) as an observable, measured at the single value 
of line-of-sight wavenumber kn — lh Mpc -1 , we find that 
our constraints on m„ for all cluster surveys improve 
by a factor of 2 if the ID power spectrum can be mea- 
sured to an accuracy, including control of systematic er- 
rors, better than AP/P ~ 1% (note this accuracy would 
allow a more modest improvement of 30% if a w = — 1 
prior was adopted). This yields an improved prediction 
of a > 2cr detection of the neutrino mass. We find that 
adding constraints from 2,000 SNe between < z < 1.7, 
whose magnitudes are measured to an accuracy of 0.15 
mag, following the Fisher matrix analysis of 32], yields 
only a modest 20% improvement on ^m„. An addi- 
tional complementary probe is tomographic cosmic shear 
statistics. Comparable bounds on m v were found for 
an LSST- like survey combined with Planck [33j| . 

In conclusion, our results suggest that detection of the 
effect of neutrino dark matter is likely to be possible. Sys- 
tematic errors such as projections and false detections in 
the WL survey may be controlled using self-calibration 
and comparisons to numerical simulations. However, we 
emphasize that additional biases will be introduced by 
selection effects in real surveys. These must be care- 
fully studied and controlled in analyzing any future sur- 
vey data. Finally, we note that at the current lower limit 
on their mass, neglecting neutrinos from the type of anal- 
ysis we outlined would result in a bias of other cosmolog- 
ical parameters. We find, for example, that w rj —0.95 
would be inferred if neutrinos are ignored in an to = -1 
universe. Hence, it is important to include neutrinos in 
any analysis that aims to derive dark energy parameters 
to ~ 1% precision. 
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